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^ i Abstract 

I We study the time evolution of the mass gap of the 0{N) non- linear sigma model in 2 + 1 

dimensions due to a time-dependent coupling in the large-N limit. We start from a thermal 
A ■ equilibrium state of the model deep in the disordered phase at initial time, let the coupling g 

. change with time towards the equilibrium critical point gc according to a fixed protocol, and study 



o 

(N 



the behavior of the mass gap under such a dynamics. Using the Schwinger-Keldysh approach, 
we derive a set of equations at large which determine the time dependent gap in terms of the 



(N 

> 
oo 

■ coupling. These equations lead to a criterion for the breakdown of adiabaticity for slow variation of 

. the coupling, which differs from the standard Landau criterion, but results in a Kibble-Zurek scaling 



law appropriate for mean field transitions with z = 2. We provide explicit numerical solutions of 
these large- equations, and demonstrate (for a protocol which has a linear ramp in the vicinity of 
^ , the equilibrium critical point gc) that there is a value of the coupling g = gc^^ > gc where the gap 

I function vanishes, possibly indicating a dynamical instability. We study the dependence of gc^^ 

on both the rate of change of the coupling and the initial temperature. We also show, by studying 
the evolution of the mass gap subsequent to a sudden change in g, that the model does not display 
thermalization within a finite time interval to and discuss the implications of this observation for 
its conjectured gravitational dual as a higher spin theory in AdS4,. 

PACS numbers: 
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I. INTRODUCTION 



The study of non-equilibrium dynamics in closed quantum system has received a lot of 
attention in recent times, largely because of its experimental relevance in cold atom physics 
l|-l8|. Typically such studies attempt to understand the properties of observables of a 
quantum system whose time evolution is governed by a time-dependent Hamiltonian, e.g. a 
time dependent coupling. Apart from their experimental relevance, this class of problems 
are theoretically interesting in their own right for several reasons. One of the reasons relates 
to their role in understanding the long-time behavior of closed integrable and non-integrable 
quantum systems. In particular, a lot of recent studies have focussed on the properties of 
long-time steady state of a closed quantum system following a quantum quench. The other 
theoretical motivation is a study of the dynamics of a quantum system during or subsequent 
to a linear or non-linear ramp of one of the parameters of the Hamiltonian. Such a dynamics 
becomes particularly interesting if the quench approaches or crosses an equilibrium critical 
point. In this case, adiabaticity will be inevitably lost close to criticality and the subsequent 
dynamics will carry universal signatures of the critical point. For slow djTiamics, a simple 
scaling hypothesis can indeed be used to show that several quantities such as density of 
excitations and excess energy scales with the rate of quench according some universal power- 
law which is determined solely by the universality class of the critical point [l| . However, in 
spite of all the advances made in recent years, there are few theoretical tools to investigate 
these issues. Remarkable exceptions in this regard are systems in 1 + 1 dimensions, where 
methods of boundary conformal field theory can be used to obtain results for correlation 
functions for an abrupt quench from a massive theory to a critical point [2^6|. However, 
(at least in higher spatial-dimensions, and for dynamics involving finite-rate ramps) unlike 
equilibrium critical phenomena, there is no general conceptual framework which explains a 
possible universal behavior of physical observables. It is thus important to accumulate as 
many "data points" as possible by examining individual models. 

The equilibrium properties of the 0{N) non-linear sigma models (NLSM) have been 
widely studied in several contexts. One such example involves a class of spin-models, namely 
Heisenberg models with antiferromagnetic ground states, which have wide relevance as ap- 
propriate models describing low-energy properties of a large class of spin-systems. A specific 
example involves double-layer antiferromagnets; it is well known that their low-energy prop- 
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erties can be described by a non-linear sigma model with the mapping rij ~ (Sij — 821) |9|, 
where n denotes the fields of the NLSM and Si(2) denotes spins on the two-layers. Another 
example of the use of such NLSM lies in description of effective low-energy field theory of 
disordered electron gas [10]. In high energy physics, nonlinear sigma models appear in many 
contexts. Examples include the low energy theory of Nambu-Goldstone modes of a spon- 
taneously broken symmetry like chiral symmetry, and worldvolume descriptions of strings 
and higher dimensional extended objects. However, in spite of its wide application in several 
branches of physics, the non-equilibrium dynamics of this model caused by a time dependent 
coupling has not been studied so far. 

In this paper, we study the non-equilibrium dynamics of the 0{N) NLSM at large 
using Keldysh method in 2 -|- 1 dimensions. We consider a time dependent coupling g{t) = 
gi + {g f — gi) tanh^(ft) (where v is the ramp rate and gf and gi are initial and final values of g) 
such that the system is initially, at t = 0, in thermal equilibrium at a temperature T in the 
disordered (paramagnetic) phase, and study the behavior of the mass gap of the model during 
and/or subsequent to this dynamics. As is well known, because of large A^ factorization, 
the N = 00 limit of such a theory is a classical theory of 0{N) invariants. Thus, a non- 
perturbative quantum problem becomes a classical problem in terms of suitable variables. 
In this work, we explore if a formulation of the problem of non-equilibrium dynamics in this 
classical theory leads to any new insight. The central results that we find from such a study 
are as follows. First, we derive a condition of breakdown of adiabaticity of such a model 
and show that it differs from standard Landau criterion used to estimate non-adiabaticity in 
many-body quantum systems. We show that this modified criteria nevertheless leads to the 
the standard Kibble-Zurek scaling. Second, we show that the gap function always vanishes 
at a time to where the instantaneous value of the coupling, g{tQ) = g^^'^ is larger than the 
equilibrium critical point at T = 0, gf^ . We use the normalization of [9J where gf^ = 1. 
We chart out the dependence of g^^^ on the ramp rate v and the initial temperature T by 
providing a numerical solution to the time-dependent gap equation of the NLSM within 
large- A^ approximation. We then compare these results with those of a related problem - 
that of quench in a Landau- Ginzburg theory. Third, we study the evolution of the mass 
gap of the model subsequent to a sudden change in g and show that the system does 
not exhibit thermalization up to a time to till which we can numerically track such an 
evolution. We point out that whereas this observation does not allow us to make a statement 
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about long-time thermalization in the model, it certainly shows that the model does not 
exhibit thermalization for t < to. We discuss the consequence of such an absence of "fast" 
thermalization (for t < to) of the gravity dual of this model. 

Large-A^ quantum quench has been indirectly studied by using the AdS/CFT correspon- 



dence 



11| to map 



these examples 
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arg e-iV limits of strongly coupled field theories to classical gravity. In 



14| (which involve N x N matrices, rather than component vectors). 



it is almost impossible to solve the field theory itself, but its dual gravity description is 
tractable. In contrast, the dual formulation of the 2+1 dimensional 0{N) vector model has 
been conjectured in [l5| to be a higher spin gauge theory in AdS^ [16] which contains an 



infinite number of massless 
papers, see in particular |17l . 



ligher spin fields. This conjecture has been explored in may 



18| 291]. In this case the field theory is tractable, and it will be 



interesting to see if this teaches us anything about quench and in particular thermalization 
time of the higher spin theory, specifically because there is an explicit dual map [17]. 

Quantum quench in the large- expansion has been studied earlier for the linear sigma 
model in [j] for infinitely fast quenches. This work does not deal with the issue of scaling 
behavior near the critical point. A similar work deals with BCS theory with an abruptly 



changing coupling 



20|. In contrast, our work in the non/mear model concentrates on the 



dependence of quench dynamics on the rate of quench. The fact that gdyn is larger than 



2l| and 



gc is similar in spirit to the phenomeon of stimulated superconductivity found in 
studied in the AdS/CFT context in 22 1. 

The plan of the rest of the paper is as follows. In Sec. HIl we use the Keldysh formalism 
to derive a set of equations which determine the time-dependent gap (which we call the 
gap function) in terms of the time dependent coupling. This is followed by Sec. IIIII where 
we derive the condition for breakdown of adiabaticity for the model as we approach the 
equihbrium critical point. Next, we discuss the determination of gdyn of this model in Sec. 
HVl In Sec. |V] we compare our results with those in a Landau-Ginsburg model. In Sec. IVH 
we discuss the time evolution of the mass gap subsequent to a sudden change of coupling. 
Finally, we conclude in Sec. IVIII 
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II. THE MODEL AND THE GAP EQUATION 



The action of the model is given by 

N 



S = cPxdt 



t)) ■ , t)) + \{x, t) [0 ■ - 1] 



(1) 



where t) is a N dimensional vector with real components. In the large N limit, N ^ oo 
with g{t) remaining 0(1). 

The field X{x, t) is a Lagrange multiplier which imposes the constraint t) ■ t) = 1. 
Redefining fields 

-)■ ^/^ = 



1 



the lagrangian density becomes, up to a total derivative 

2 V 9[t)J 9{t) 



where 



3 (i 

2 \g 



t) = A(x, t) + a(t) 



(2) 



(3) 



(4) 



The last term in ([3]) is field independent and can be therefore ignored. 

The partition function can be expressed in the Schwinger-Keldysh formalism as 



(5) 



where we have doubled all the fields as usual, and S is the action for the lagrangian in (j3]). 
As is well known, the representation (|5]) is schematic 10|- one has to pay attention to the 
end point of the time contour. However these "boundary" terms do not affect the saddle 
point equation, though these are important for evaluation of the partition function by the 
saddle point solution. 

One can now integrate out the fields ip± leading to the effective action for T,±, 



Seff = NTr log(D-i) -N J £xdt-^[J:+ - E_] 



(6) 



where D is the propagator matrix whose inverse is 



-92 + s_ 



(7) 



The large-N saddle point equations therefore become 



^ TrD++, J_ = Tr(-D__) (8) 



9{t) 9{t) 
Therefore, the saddle has S+ = S_ = and the equation becomes 

1 



, cPx <x,t\D\x,t>3 (9) 

g{t) J '11'/^ y J 

where we are considering the problem at a temperature T = Note that the equality 

of S+ and S_ is a feature of the strict N = oo limit. Fluctuations around the saddle point 
will destroy this equality. 

The coincident Green's function in may be obtained by considering a Heisenberg 
picture field which satisfies the homogeneous equation 

[-d^+d^-mMx,t) = o. (10) 

Using a mode decomposition 

/(Pk 
j^[akXk{t) + alxlit)] (11) 

the equal time Green's function < x(x, t)x{x', t) >i3 is the two point function in the thermal 
state, i.e. in the state 

where 



The solution xa;(^) may be written in the form 



where flk{t) satisfies the equation 

/ " \ ^ 

^' +nl = P + E{t) (15) 



2 

2 fifc 4 I 



Then the gap equation becomes 



(?(t) 7 (27r)2 2fifc(t) 



coth(^) (16) 
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Note that the exponential factor in (fT^ canceled in the expression for the coincident time 
Green's function. 

The equation (ITB|) has to be solved for Qk(t) for a given g(t) and substitution of the 
solution in f[T^ gives the gap function 

Before ending this section, we note that for a time independent coupling qq, our formalism 
reproduces the well-known equihbrium solution {q]. In this case, equation (fT5|) shows that 
fifc = A;^ + S i.e. the gap is now independent of time. The integral over k on the right hand 



side yields the result 
Here gdf^) is given by 



2tt 



1 

9o 



(3 



log 



2 sinh( 



/3 



log 



2 sinh( 



and A is the momentum space UV cutoff. The equation ( fT71) can be solved for So 



2 2 ^ 



(17) 



(18) 



(19) 



For any non-zero temperature this equation has a solution. At exactly zero temperature the 
gap equation becomes 

/^o = (20) 

so that a real solution exists only for > g^. For g^ < g^ the 0{N) symmetry is spon- 
taneously broken and the theory is massless. The large-N solution presented above is not 
valid in this phase, though a valid solution in this phase is well known jol. The point go = gc 
is then a critical point which separates the ordered and the massive phase. At non-zero 
temperatures, the critical point becomes a cross-over which separates regions of the phase 
diagram which are qualitatively similar to the ordered and the disordered phases. The 
location of the cross-over point is given by ~ 1 which implies T for small T. 



III. QUANTUM QUENCH : BREAKDOWN OF ADIABATICITY 

When the coupling g{t) varies slowly compared to the mass scale set by the coupling 
itself, one expects that the gap function S(t) evolves adiabatically. Adiabaticity should 
break down when the gap function becomes small, e.g. near the zero temperature critical 
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point. Let us define 

h{t) = 27r(^- - 4t) (21) 

To investigate adiabaticity, we first need to find an expansion for VLk{t) in terms of a 

by solving (JT5l) in a derivative expansion. Tliis is easily done, and the lowest order result is 



^^fc(t) VP + 8(A;2 + S(t))5/2 32(P + S(t))V2 

We then need to substitute this in (fT6|) . 

For zero temperature, it is possible to perform the necessary integrals and the lowest 



order result is 



IS IS 



2 



ftW = v^-5j^ + 55^ + --- (23) 



Inverting this (again in a derivative expansion) we get 



Therefore adiabaticity breaks down when 



2| _ 1,2 



h' (25) 



'h 

We will be interested in generic profiles of g{t) for which h{t) ~ vt near t = 0, e.g. h{t) = 
— atanh(ft) for some dimensional parameter a. For such a profile adiabaticity breaks at a 
time 

t. ~ V-'/' (26) 



which is the usual Kibble-Zurek scaling for linear quenches [30| . We note that a naive use of 
Landau criterion in this case would lead to a different scaling law. According to the Landau 
criterion, the adiabaticity breaks when the gap A(t) = -\/S obeys A ~ A^. This leads to 



A(t) ~ ^/hit) ~ (vty/^ ^ ~ (27) 

which is not the usual Kibble-Zurek scaling. 

These result trivially extends to nonlinear quenches. We note that the breakdown of 
adiabaticity can be also investigated analytically for low temperatures and the results are 
qualitatively similar up to exponentially small corrections. 
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IV. THE DYNAMICAL INSTABILITY 



In this section, we determine the value of the couphng where the gap function first 
becomes zero. This possibly signals a dynamical instability. The protocol for h{t) that we 
follow for studying this phenomenon is the following: we start with a fixed g = Qi inside the 
disordered phase and with an equilibrium temperature T and decrease g to gf < Qi a.t the 
end of the evolution with a speed v. This protocol is realized using 



gii) = gi + {gf - gi) tanh^(f t) 



(28) 



In what follows, we track the time evolution of the gap and focus on finding the largest 
value of gj for which the minimum value of the mass gap, Smin, reaches zero at some point 
during the evolution. This value yields g"^^^. Since the numerical solution of the gap equation 
becomes difficult when Smin — 0, we extract the position of the dynamical critical point by 
extrapolation of as a function oi gj. 
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FIG. 1: Rate dependence of gc^° — 1 for T = 2 (red circle) and Qc^^ for T = 1.75 (blue circle). The 
shift in the value of gt^^ for T = 2 is carried out to enhance clarity. 

The numerical procedure we adopt is the following. First, we provide initial guess values 
of E(t) for a discrete set of points t = ti and numerically solve Eq. [15] to obtain Q{ki,ti) 
for an array of discrete fcj and ti. From these values, using interpolation, we compute the k 
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integral appearing in the right side of Eq. [16] and obtain a set of trial values (?triai(^j)- We 
then minimize the function [g{ti) — gtviaiiti)]"^ self-consistently by varying We choose 

the size of the set of points tt and ki such that variation of size of the sets do not change 
numerical values of 

The results of this procedure is summarized in Figs. [T] and [2J We find that for all rates 
V > 0.25 and temperatures T > 0.05 that we have studied, the gap function first touches 
zero at g = g'^^^ > gc- Note that in our convention gc = 1. The fact that g^yn is larger than 
gc is similar in spirit to the phenomeon of stimulated superconductivity found in 2l| and 
studied in the AdS/CFT context in 22 1. 



We find from Fig. [T] that the value of g^^"^ increases with increasing v, reaches a maximum 
around a critical rate v* which depends on the starting equilibrium temperature, and then 
decreases as v is further increased. We find v* reduces with decreasing temperature and the 
hump flattens. 

The existence of v* can be qualitatively understood as follows. Since very slow dynamics 
in the disordered regime is expected to be adiabatic, we expect g'^^^ to approach gc for small 
V. On increasing v, g^^^ increases and deviates from gc- This continues till a rate v = v* after 
which the system does not have enough time to respond to the drive leading to a decrease 
in (yf^^° with increasing v. Note that here we have restricted our numerics to values of v so 
that the system reaches gj around t = Iq till which we track the dynamics. For faster v and 
fixed to; the system will eventually enter the quench regime where g^^^ will be determined 
by the evolution of with g = gf subsequent to the quench till t = to. We do not address 
this regime in this section. 

The temperature dependence of g^^"^ for different v is shown in Fig. [2l Here we find that 
there is a crossover regime around T = 1 where the behavior of g^^^ changes from the high 
temperature region where it shows appreciable variation with f to a low temperature regime 
where it becomes virtually independent of v. Note that the saturation of g^^"^ for small v is 
consistent with the shifting of v* to lower values with decreasing temperature. 

We end this section by noting that the dynamical instability of the time-dependent mass 
gap studied above need not correspond to a phase transition since the nature of the correla- 
tors of theory at the point where the time dependent mass gap vanish need not be long-range 
due to memory effects incorporated in Qk{t)- It would be interesting to study such correla- 
tors near the instability; presently the slow convergence of the numerical solution near the 
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FIG. 2: Temperature dependence of gc^" for v = 2 (black circle), 1.5 (red square), 1.0 (blue 
triangle), 0.75 (violet hexagon), and 0.5 (green inverted triangle). 

instability prevents us from carrying out such a detailed study. 



V. A DIGRESSION : LANDAU-GINZBURG DYNAMICS 

It is instructive to compare the above behavior with the non-equilibrium dynamics of a 
toy Landau- Ginzburg (LG) model. Consider a LG model with z = 1 m the presence of a 
time dependent source J(t), with the mass being tuned to zero. The equation for the order 
parameter is given by 

The equilibrium critical point is at J = 0. 

The motivation behind considering this toy model is the following. Large N theories are 
classical in the leading order, and one might imagine that there is a classical equation of 
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motion which describes this hmit, typically in one higher dimension. Models of N x N 
matrices are often of this type, e.g. matrix quantum mechanics whose large-N limit is 



described by the classical equations of 1 + 1 dimensional string theory [23|] - the string field 
in this case is in fact a single massless scalar which can be identified with a suitable collective 
variable AdS/CFT dualities are also of this kind : the large-N classical theories are 
generically string theories in one higher dimension and contain an infinite number of fields. 
Nevertheless, usually in the strong coupling limit, only a few of these fields are massless 
(which is the statement that in the field theory only a few operators have dimensions of 
order one rather than of order A^). In that case the dual theory is described by a finite 
number of classical equations of motion, viz. Einstein equations and equations of motion of 
a few other fields. Our toy model is supposed to be carricature of such large-N limits. In 
fact, the dynamics near phase transitions in AdS/CFT models can be accurately captured 
by Landau-Ginsburg theories like (129|1 jl^ . 

The large-N theory of the 0{N) model as formulated above does not appear to be of this 
type - the classical saddle point equation is not a differential equation of the LG type, but 
an integral equation. Nevertheless, for a slowly varying coupling, the gap function will obey 
an inhomogeneous differential equation at large N, where h{t) appears as a source term. In 
fact to lowest order the equation is given by (1231) . Of course this kind of equation fails when 
the gap function becomes small. Our purpose is to compare the behavior of the 0{N) model 
and a LG type model. 

Consider a time dependent source in (l29l) of the form 

J{t) = Ji + {Jf - Ji) tanh^ivt) (30) 

pretty much like g{t) chosen for the 0{N) model. We now solve the equation for a fixed Ji, 
starting at t = with adiabatic initial conditions for different values of J/, and determine the 
value of Jf = Jj for which the order parameter ^ first touches zero. A typical time evolution 
of the order parameter is shown in Figure ([3]) for vanishing friction and in Figure (jlj) in the 
presence of some friction. 

Figure ([5]) shows the behavior of Jdyn as a function of the ramp speed for vanishing 
friction. Figure ([6]) is the same plot in the presence of friction. For very small v, Jdyn is very 
small and increasing with v since the time evolution is expected to be adiabatic. For large 
V we have a rapid quench. The system does not have enough time to react to the change 
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(c=0, J_i=4, v= 0.24, J_f = 0.00215) 

Order Parani eter 




FIG. 3: A typical behavior of the order parameter as a fmiction of time, with the source function 
given by ()30p and vanishing friction. The evolution starts from t = 0. The purple line is the 
adiabatic value of the order parameter. The value of Jf chosen in this plot corresponds to the 
dynamical critical point. 

ic=0.1, J_i=4, v= 0.66, J_f = 0.02851 
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FIG. 4: A typical behavior of the order parameter as a function of time, with the source function 
given by (jSOp and some friction. The evolution starts from t = 0. The purple line is the adiabatic 
value of the order parameter. The value of Jf chosen in this plot corresponds to the dynamical 
critical point. 

and remains in the initial state for quite some time, after which it starts oscillating. Jdyn 
saturates to a constant value. 

For vanishing friction, the saturation value may be understood as follows. For large v, 
the quench is rapid at t = 0. We may then approximate J(t) by J(0) = Ji and J(t) = Jf 
for t > . Then a first integral of the equation of motion fl2^ for t > is given by 
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(J=J_i+(J_f-J_i)*(tanh(vt)'^2;J_i=4,c=0) 

log(L(dyn)) 



FIG. 5: A log-log plot of Jdyn versus v with no friction. The saturation value is 4 Jj as expected. 



{J=J_i+(J_f-J_i)*(tanh(vt)'^2;J_i=4,c=0.1) 

log(]_(dyn)) 



log(v) 



FIG. 6: A log-log plot of Jdyn versus v with friction. 



Since at t = the initial conditions are adiabatic ^(0) = [-Ji^^^ and ^(0) = 0, so that 
E = jjf^^ — Jfj]^^. Now, at the time when the order parameter first touches zero, we must 
have ^ = and ^ = 0. For this to happen we must have = 0, i.e. Jj = \Ji. Indeed, we 
have checked that the saturation value of Jdyn is indeed given by ^ Jj. 

The behavior of Jdyn as a function of f , however, displays non-monotonic behavior for 
intermediate v, with multiple well-defined humps. In the presence of friction we could not 
find multiple humps, but one hump always remains. We do not understand the reason behind 
this behavior. However, it is interesting that there is a similar behavior in the nonlinear 
sigma model, even though the latter differs significantly from our toy model when the gap 
vanishes at some time. 



14 



It is also interesting to consider a profile of J{t) given by 

J{t) = {a~b) + btanh^{vt) (32) 

and consider the time evolution starting at a large negative value of t proceeding to a large 
positive value of t. For a given a the minimum value of J{t) is (a — 6) at t = and this come 
close to the equilibrium critical point as b comes closer and closer to a. For a fixed value 
of a we then solve (1291) for various values of b and find the b for which the order parameter 
touches zero at some value of J = Jo{v). A typical behavior of the order parameter as a 
function of time is shown in Figure ([7]) 



I Order Parameter, v=0.5, a =5, b=4.86| 




FIG. 7: The order parameter as a function of time for a time dependence given by (|32p . The purple 
line is the adiabatic solution 

It may be verified that so long as the initial value of the source, a, is large enough the 
location of Jo does not depend on a. We then study the behavior of Jo as a function of v. 
The result is shown in Figure® in a log- log plot. 

For small values of the rate v, the behavior is Jo ~ f For large values of v the Jo 
saturates. 

The small v behavior is consistent with the scaling properties of the solution to (1291) . 
When V is small, the value of (a — b) for which the order parameter hits a zero is small as 
well, so that one may approximate J(t) ~ {vt)"^- In this case, the solution to this equation 
may be easily seen to be of the form 

^{t,v) = v^/''(j){tv^/^l) (33) 

which shows that at the typical time scale t ~ f ^^/^ the value of the source is J ~ v^^^ = v^''^. 
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{ J=(a -b) +b(tanh (vt)"2;a =5) 

log(LO) 
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FIG. 8: log{Jdyn) as a function of logv. For small v one gets Jd- 
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The presence of humps in the behavior of Jdyn is intruigingly similar to similar non- 
monotonicity found in the 0{N) model. However, other aspects of these results appear 
to be quite different from our 0{N) results. In particular, while Jdyn or Jq iii the LG 
theory saturates at large speeds, the g'^^"' appears to first rise and then come closer to the 
equilibrium critical coupling Qc = 1. It is conceivable that for much larger speeds, g^^"' rises 
again and saturates. However, numerical convergence becomes difficult at high ramp speeds. 

It is also conceivable that this difference reflects a fundamental difference between the 
large-N classical theory of 0{N) as formulated above and a LG type theory. As discussed 
above the gap function is not related to the coupling in a local fashion in time. The gap 
equation ( !T6|) relates g(t) with flk{t) which depends on the entire function through the 
differential equation ( ITSl) . It is only in the adiabatic approximation that the gap function 
satisfies a differential equation with a source given by the coupling - and adiabaticity of 
course fails when the gap function vanishes. 

There is, however, a formulation of the large-N limit where the saddle point equations 
are in fact usual differential equations. In fact the complete set of 0{N) invariants which 
reproduce singlet sector correlators is given by bilocal fields 

a{x,t-y,t) = j^^\x,t)^\y,t) (34) 



The collective field theory of these bilocal fields 



27| can be written standard field 



theory of an infinite number of massless hig. 
spectrum as Vasiliev's higher spin theory 



ler spin fields [17[ which has precisely the same 



16| . thus providing an explicit construction of 



AdS/CFT duality in this case. In this formulation, the gap function is simply one of these 

16 



fields, and the effective action we use is the result of integrating out everything else. However, 
unlike other AdS / CFT examples, there is no scale which separates the fields which are being 
integrated out - so what we get is a nonlocal field theory of the gap function alone. It is 
tempting to speculate that this basic difference is somehow responsible for the difference in 
the dynamics, though we do not have a precise understanding of this issue. 



VI. TIME EVOLUTION AFTER A SUDDEN QUENCH 

In this section, we discuss the time evolution of the mass gap subsequent to a 
sudden impulse imparted to the system. The impulse is imparted by changing g{t) as 
g{t) = [gi — gd) + gdiaj^[v{t — ti)] with gi = 4, (7^ = 1, ti = 1 and v = 20. Note that 
with this choice, the system starts its evolution with g = g^ sX t = 0. Near t = ti, the 
coupling g changes to gf = gi — gd and back to g^. The change takes place within a time 
window of r ~ 1/f around ti and thus appear as an instantaneous impulse for large v. The 
plot of the subsequent evolution of S is plotted as a function of time in Fig. M for several 
initial temperatures. We note that for all temperatures, the system does not show any sign 
of thermalization in the sense that does not approach any constant steady-state values 
till the time to = 10 that we track it's evolution numerically. 

We believe that this is a manifestation of the lack of thermalization in 0{N) vector 
models to leading order in We note that from our numerical result, we can not 

rule out thermalization of the system at longer times; however, the model certainly do not 
thermalize for t < to- Such a lack of thermalization in response to an impulse in this model 
appears from the lack of quasiparticle scattering at 0{1/N). Such scattering in the present 
models appears in 0(1/A^^) and ultimately leads to thermalization. This is in contrast with 

n 

large-N models of matrices, where thermalization is expected to occur [25| at the leading 



order 



31| This is manifest for the class of such large-N models which have gravity duals. 



12| . In such models, thermalization 



where thermalization is seen as formation of black holes 
is almost instantaneous for local operators. 

This is also consistent with the fact that the gravity dual of the 0{N) model is a higher 
spin gauge theory rather than standard Einstein gravity. In usual AdS /CFT duals of models 
of large-N matrices (e.g. gauge theories), thermalization is signalled by black hole formation. 
However, a study of the finite temperature properties of the singlet sector of the higher spin 
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model shows that there is no large-N transition at order one temperatures [2u\ . This possibly 
implies the absence of thermodynamically stable large black holes with order one Hawking 
temperatures in this higher spin theory. 



2.0 n 




FIG. 9: The evolution of the mass gap after a sudden quench at t = 1 as a function of time till 
t = to = 10 for temperatures T = 2 (blue dotted line), T = 1 (red dashed line) and T = 0.1 (black 
solid line). 



VII. CONCLUSIONS 



In conclusion, we have studied the non-equilibrium dynamics of large vector model in 
it's paramagnetic phase using Keldysh formalism. Although in this work we have primarily 
concentrated on the behavior of the model in 2 + 1 dimension, the formalism developed 
may also be used for one and three spatial dimensions. We have shown that the condition 
of breakdown of adiabaticity in this model due to a time dependent coupling g(t) differs 
from the usual Landau criterion, but the analysis reproduces the standard Kibble-Zurek 
law for usual mean field transitions with z = 2. We found that during the time evolution, 
the instantaneous value of the gap reaches zero for a coupling gdyn which is larger than the 
coupling at the equilibrium critical point and studied the dependence of gdyn on the quench 
rate and compared the results with those in a Landau-Ginsburg dynamics. Finally, we have 
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studied the response of the model to a sudden perturbation, shown that the model does not 
exhibit short-time thermalization, and have discussed the consequence of this phenomenon 
for it's gravity dual. Our work can in principle be extended to studying non-equilibrium 
dynamics of the symmetry broken phase of the model, to dynamics of the model in 3 + 1 
dimensions, and also to its periodic dynamics. We leave these issues for future studies. 
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